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I. INTRODUCTION 


The Schrodinger equation is the starting point to taking into account the quantum me¬ 
chanical effects in physical systems. The Klein-Gordon (KG) and Dirac equations have to be 
investigated if one also would study the relativistic effects on them. It is well known that the 
potential forms for which KG and Dirac equations can be solved exactly are very restricted. 
So, the approximate solutions of the above equations for equal and/or unequal scalar and 
vector potentials become important to achieve some informations about physical system. 
The potential helds having exponential shapes are one of the important part having exact 
or approximate solutions for KG and Dirac equations, and have received great attention in 
the last few decades [1-15]. 

The case where scalar and vector potentials are nearly equal in magnitude with oppo¬ 
site sign, S'(r) ~ —U(r), is called the pseudospin symmetry, and the Dirac equation has 
pseudospin symmetric solutions in that case [16, 17]. If scalar and vector potentials sat¬ 
isfy the condition S{r) ~ V{r) then the Dirac Hamiltonian has the spin symmetry [18]. 
The pseudospin symmetry is an exact symmetry for Dirac Hamiltonian under the condi¬ 
tion d{S{r) + V{r))/dr = 0 while the spin symmetry is an exact one under the condition 
d{—S{r) -|- V{r))/dr = 0 [19, 20]. We present here the approximate analytical solutions of 
Dirac equation for the case where S{r) ~ —V (r) corresponding to pseudospin symmetry for a 
potential held which could be written in a exponential form. This hyperbolic-type potential 
has been studied for the non-relativistic case [21, 22], and the effects of position-dependent 
mass on spectrum have also been computed [23]. The explicit form of the potential (Fig. 1) 
is 


V(r) = -Vc 


sinh^{r / 2d) 


( 1 ) 


^ cosh^{r /2d) ’ 

where the potential parameters Vq and d correspond to depth, and width of potential, 
respectively. Becasue of the p and q the above expression dehnes a class of potentials such 
as g = —2, 0, 2,4, 6 and p = —2, 0,... g. The case (g,p) = (2, 0) describes the Poschl-Teller 
potential. For the potentials with g = —2, 0, 2, the Schrodinger equation can be converted 
into the type of hypergeometric equations [21]. The Schrodinger equation for the case where 
g = 4, 6 could be written in the form of a confluent Heun equation [24]. In the present 
work, we achieve also a confluent Heun differential equation [25, 26] by converting the Dirac 
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equation for the case where {q,p) = (6,4) by using suitable transformations on coordinate 
variable r. 


The plan of this paper is as follows. In Section 2, we give the basic equations briefly which 
also contain some mathematical properties of the solutions of confluent hypergeometric- 
type equation. In Section 3, we handle a transcendental equation which could be solved 
numerically for the bound states. This equation is written in terms of the wave functions, 
and we provide the zeros of the real and imaginary parts of transcendental equation which 
correspond to the points on ” E-axis” where the two curves cross. We hnally give our 
conclusions. 


II. THE OUTLINE OF DIRAC EQUATION 


The Dirac equation for a particle with mass p' is written as [27] 
[a.p + /3[/i' + S{r)] + V{r) - E] T(r) = 0 , 


( 2 ) 


where E is the energy, p is the linear momentum, a and /3 are 4x4 matrices have following 
form, respectively, [27] 


d = 








Here, a is Pauli matrix, and I corresponds to 2 x 2 unit matrix. The spherically symmetric 
Dirac spinor is written in terms of the quantum numbers n (usual quantum number), k (spin- 
orbit quantum number), i (orbital angular momentum quantum number), m (x-component 
of angular momentum quantum number), and has the following form 


if Eir)Yl^{e,<j>)\ 


( 3 ) 


where ^ = ^ + 1, and k = ± (j -|- 1/2) with j = £ ± 1/2 (total angular momentum quantum 
number). The functions Y^^iO^cf)) and Yj^{6,<p) are spherical harmonics, and E{r)/r, and 
G{r)/r are radial part of the Dirac spinor. By inserting Eq. (3) into Eq. (2) we have two 
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coupled differential equations 


^ ^ I F(r) -[^I' + E- V{r) + 5(r)]G'(r) = 0 , 


d K 
dr r 


G{r) -E + V{r) + ^(r)]F(r) = 0 , 


(4a) 

(4b) 


By using Eq. (4a) for F(r), and substituting it into Eq. (4b), we obtain two independent, 
second-order differential equations as 

{ ^ + (f^'-E + C)[V(r) - S(r)] } G(r) = + C(^' + E)]G(r) , 

(5a) 

{ ^ -(t^' + E- G')\V(r) + S(r)] } E(r) = [,,« - £= - - i!)]F(r), 

(5b) 


where the constants C = E(r) -|- S{r) (i.e., dC/dr = 0), and C = E(r) — S{r) (i.e., 
dC/dr = 0), respectively. In the next section, we deal with the solutions of Eq. (5a) for the 
potential under consideration. 


III. HYPERBOLIC-TYPE POTENTIAL 


Inserting the potential with (6,4) into Eq. (5a), taking into account V{r-) + S{r) = C, 
and using the approximation instead of the centrifugal term 

1 ^ 1 

Asinh‘^{r/2d) ’ 


gives us 


dx"^ 


Ad^ 


( k{k — 1 ) 

\Asinh‘^{x) 


MEn 


sinE^[x) 


G{x) = 0, 


( 6 ) 


cosh^{x) 

where x = r/2d, M = fi' — E + G, and e = — E"^ + G+ E). 

By making a new transformation y = 1/cosK^{x), writing the wave function G{y) = 
e^yf{y), and with the help of following abbreviations 


ai = —^ d?‘K{K — 1) — d^MVo — d/^e , 
02 ~ cq d MVq , 
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(7a) 

(7b) 



we obtain a differential equation for f{y) as 




#2 


y ^-yjdy 


d 

—I" 


A -j- fli 

y 


+ 


0,2 — A /2 
1 -y 


d^e d'^K(K — l)/4 
y^ (l-|/)2 


0 , 

( 8 ) 


with = —d^MVo- 

Finally, the transformation such as f{y) = y^'{l — yY'H{y) gives us a second-order 
differential equation 




1 + 2 p' _ 2q' + 1/2 

y 1 -d 


A 

dy 
p \ 1 


2Ap’ + A + oi- 2p'q' - O'' - 77 - + 

2 ) y 


02 


j - 


2p'q' 


/ P 

q - 

2 



0 , 

(9) 


where the constant p' and q' in this equation should satisfy the equalities p'{p' — l)+p' — d‘^e = 
0, and q'{q' — l) + q'/2, — {l/4:)d‘^K,{K, — 1) = 0, respectively, for getting a confluent Heun-type 
equation [25, 26]. If we write Eq. (9) as 


d^ 

dy^ 


+ 2A + 


l + Bi 1 + B2\ d B^y + B. 


y 


+ 


+ 


y-l J dy y{y-l) 


>H{y) = Q, 


( 10 ) 


where 


Bi = 2p' ; B 2 = 2q' - ( 1 / 2 ) ; B^ = 2A{p' + q') + 01-02 + 3A/2 

Bi = {1 + 2p'){q' — A) — oi -|- pY 2 , ( 11 ) 


and compare with the form [26] 


dp"^ 


+ 


d+17+l 

Ol H-h 

p 



ff(p)=0, 


( 12 ) 


we write the solution of Eq. (10) as Hiij) = HeunC{a, (3,7 , d, 7 , y). Eq. (10) has two regular 
singularities at y = 0, and 1, and an irregular singularity at y —)■ cx). So, the solution is 
given by [26] 

CXD 

H{y) = HeunC{a, /3, 7 , S, p,y) = ^ dny'^ , (13) 

i=l 

where the coefficients d„ satisfy the relation A'^dn = B'^dn-i + C'^dn -2 with the conditions 
do = 1 and d_i = 1 with initial conditions d_i = 0 and do = 1. The coefficients in this 
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summation are 


^n = l + /^/^, (14a) 

^n = l + (/^ + 7 - «- 1)/^ +{v- /^)/2 + (7 - a)(/3 + l)/n2 , (14b) 

= q;((5/q; + (/9 + 7)/2 + n - l)/n^ . (14c) 

where A'^ ^ 1, ^ 1 and —)■ 0 if n —)■ cxd. It is assumed that the normalization of 

HeunC{a, (3, 7 , S, rj, y) = 1 . The parameters introduced in last equation satisfy the relations 

/i= (l/ 2 )(a-/ 3-7 + / 3 (a- 7 ))- 7 , (15a) 

z/= (1/2 )(q; ++ 7 + 7 ( 0 ; +/9)) + (5 + 7 . (15b) 

In the present work, we have the set of equations about parameters 


a = 2A ] 13 = Bi ;7 = i ?2 ;y = --B 4 = B 3 + B 4 , 

7] = yl(l + 2Bi) + B 4 — 1 (i?i + B 2 + B 1 B 2 ) S = -B 3 — A(2 + Bi + B 2 ), 
We need two conditions as follows [24-26] 

+ 


S = -a\^N + l + 
Atv+i = 0, 


(16) 


(17a) 

(17b) 


to obtain a confluent Heun polynomial of degree iV > 0 for y (see Ref. [21], for details). 
Aat+i in Eq. (17b) is a (TV-1-1) x (iV-|-l) determinant. In the light of the above requirements, 
we write two linearly independent solutions [28] 

~ exp[Asech^{xy\sech^^ {x)tanh^‘^ {x)HeunC{a, (3, 7 , S, 7 , sech^{x)) , (18a) 

G‘'‘^\x) ~ exp[Asech?{x)]sech?‘^'^'~^^\x)tanh?'^' {x)HeunC{a, I3,'y,5,rj,sech?{x)) , (18b) 


The confluent Heun functions HeunC{a, (3,^,6, ri,y) are convergent within \y\ < 1, and 
asymptotic behaviour at y = 1 is not known [ 21 , 23-28]. So, by expanding the solution about 
the other singular point we could obtain an analytical continuation for the Heun function. 
We hnd a transcendental function corresponding also the Wronskian written in terms of 
wave functions [25, 28] which could be solved numerically for relating two solutions. The 
other set of solutions can be obtained by defining a new variable z = 1 — y, m this instance 
Eq. (10) becomes 


dz'^ 



l + 5i 
z-1 


1 + -B2 \ d ^ / B^ + B4 ^ 54 \ 
X J dz \ X X — 1/ 


>H{z) = 0, 


(19) 
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with two linearly independent solutions 


G^^\x) ~ explAsech"^{x)]sech‘^^ {x)tanh?'^ {x)HeunC{—a,'y, —S,ri + S,tcinh‘^{x)), (20a) 
G^‘^\x) ~ exp[Asech‘^{x)]sech‘^^^'~^^\x)tanh‘^'^'{x)HeunG{—a, —7,/9, —S,p + 6 ,tci-nh'^{x)), 

(20b) 

Within the quantum mechanics, we should construct mathematical equalities for the wave 
functions satisfying the continuity conditions which will be discussed in next section. 


IV. BOUND STATES 


We could require that two solutions, ^*^^^( 0 ;) and G^^\x), are sufficient to obtain the 
energy spectrum. Quantum mechanical restrictions say us that the wave functions and its 
derivatives must be continuous such as 

and 

dG^^\x) _dG^x)^ 

dx dx 

which could be combined as a determinant 

G«(a;i) G(3)(a;i) 

dGW(x) I 

dx hi dx hi 

where providing a point (xi) in analytical domains of both Gh)(a;) and G^^^(x). Eq. (21) 
gives us a transcendental equation, and we write for nontrivial solutions 

T(E) = (^hlM) I _ !-.='>• (22) 

\ dx J \ dx J 

Since the function J^{E) is complex, the independent equations Re[E{E)] = 0 and 
Im[E{E)] = 0 should be solved to obtain real solutions (bound states energies). These 
two equations can be solved numerically because they are transcendental equations and the 
energy eigenvalues are the points which correspond to the one on ” E-axis” where plots of 
Re[E{E)] = 0 and Im[E{E)] = 0 are cross. In Fig. (2), we show the dependence of E{E) 
on E with Xi = 1 for the potential parameters Vq = 200, and d = 0.5. The trigonometric po¬ 
tential has three states at points E = —6.60, —11.83 and —16.13 for p,' = 2.0, while the ones 
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at points E = —5.99, —11.26 and —15.63 for p' = 0.5, and also at points E = —5.88, —11.12 
and —15.53 for p' = 0.01, respectively. The values of bound state energies shift to left on 
” E-axis” while the mass increases which means that the energy levels appear in more deeper 
of potential-well. 


V. CONCLUSION 

The hyperbolic-type potential has been solved for a Dirac particle under the case where 
scalar and vector potentials are equal in magnitude with opposite in sign. The solutions 
have been computed in terms of confluent Heun functions. The bound state energies have 
been found numerically by obtaining the zeros of transcendental function E{E) on ” E-axis”. 
It has been given two plots of E{E) for different mass values to see the bound state energies 
clearly. 
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FIG. 1: The hyperbolic-type potential for Vq = 200, d = 0.5. 




(a) Intersections of transcendental functions for (b) Intersections of transcendental functions for 
= 2. /i' = 0.5. 



(c) Intersections of transcendental functions for 

Ai' = 0.01. 

FIG. 2; The energy spectrum is derived by solving numerically the transcendental 
equations Re[F{E)] = 0 and Im[F{E)] = 0. In the figure, we show Re[E{E)] (solid line) 
and Im[E{Ey\ (dashed line). The energy eigenvalues are given by the points on horizontal 
line where two curves cross. Parameters: /t = 1, G = 0, h = c = 1. 
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